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This paper attempts to relate turbulence theories to radio measurements on vhf 
ionospheric forward scatter circuits. To this end, the single scattering description of the 
electromagnetic response of electron density irregularities and the corresponding transmis- 
sion expression are evaluated. Statistical distributions of signal levels are found to agree 
with a scattering model. The several theories for turbulent mixing of the electron density 
are then summarized. A turbulence mixing model is compared favorably with experimental 
data on absolute signal levels and their diurnal and seasonal variations. Scattering heights 
in the ionosphere responsible for the signals are also consistent with these theories. Fre- 
quency and distance dependence scaling laws are compared briefly with the data. The 
scatter signal behavior (luring sudden ionospheric disturbances is also explained. No 
attempt is made to compare meteoric and turbulence scatter contributions to the measured 
quantities in this paper. Approximately thirty references are given. 



1. Introduction 

The chance observation of u inexplicably high field 
strengths beyond the optical horizon at microwave 
frequencies [l] 2 gave the first hint of "scatter propa- 
gation." Subsequent experiments showed that the 
received signal was weak but significantly greater 
than that predicted by round earth diffraction 
theory [2]. The field strength decreased slowly with 
distance and was observed to fade about its mean 
value several times a second. These observations 
suggested a statistical theory of scattering from 
turbulent fluctuations of the tropospheric refractive 
index [3, 4]; a theory which lias been developed and 
fitted to increasingly refined experiments [5]. 

It was inferred from the fading of ionospherically 
reflected skywave signals that similar refractive 
irregularities must also exist in the lower region (E 
and D layers) of the ionosphere. This evidence plus 
the reality of tropospheric scatter suggested that a 
similar propagation mode might be sustained at 
frequencies above the maximum usable frequency 
(MUF) by scattering from turbulent concentrations 
of free electrons in the ionosphere. Such a mech- 
anism was immediately observed at vhf for ranges 
extending to 1,500 km [6], thus providing an im- 
portant new mode of radio propagation. 

The phenomenon of ionospheric scatter propaga- 
tion has also evinced a growing scientific interest. 
The basic theoretical problem combines electro- 
magnetic theory, probabilistic descriptions of ran- 
dom signals, and turbulence theory in an essential 
mixture. On the experimental side, radio waves 
represent an accurate, reliable tool for studying 
turbulent conditions at high Reynolds numbers in 
the ionosphere. 

i Contribution from Space Technology Laboratories, Inc., Los Angeles 45, Calif. 
\n abbreviated version of this paper was published in J. Geophys. Research 
4, 2230 (Dec. 1959). 

2 Figures in brackets indicate the literature references at the end of this paper. 



This paper is presented as a linkage between the 
turbulence theories and the radio measurements. 
It attempts to point up the important assumptions 
and questionable approximations in the present 
transmission expressions by way of emphasizing the 
accuracy of this linkage. At the same time, brief 
accounts of the experimental data and theoretical 
models are given in the interest of acquainting the 
radio physicists and turbulence theoreticians with 
(he implications of the basic results in the opposite 
field. The present paper is preliminary in that it 
anticipates some results of a thoroughgoing com- 
parison of theory and experiment for ionospheric 
scatter propagation, which is to be prepared shortly. 
This paper does not consider the contribution of 
small meteors to the background scatter signal. 

2. Statistical Behavior of the Signal 

Random fading of the received signal is one of the 1 
principal reasons for identifying ionospheric scatter 
propagation with electromagnetic scattering from 
turbulent irregularities. A typical time history of 
the amplitude R of the complex voltage x+iy induced 
in the receiver by the incident electromagnetic field 
is shown in figure 1. Statistical analysis of actual 
records show that the probability that the amplitude 
exceeds a given level is approximately described by a 
Rayleigh distribution [7, 8]. This suggests that the 
process R(t) is the vector amplitude of two orthogonal 
(phase quadrature) Gaussian random processes. On 
the other hand, if the received signal were produced 
by a large number of independent scattering contri- 
butions, one would expect the signal components to 
be distributed in a Gaussian manner by the Central 
Limit Theorem [9]. This line of reasoning indicates 
that scattering theory provides a reasonable (though 
not unique) explanation for the statistical distribu- 
tion of the experimental data. 
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Figure 1. Typical amplitude versus time record of ionospheric 
scatter propagation signal amplitude together with correspond- 
ing time — displaced vector voltage diagrams which are assumed 
to describe the process. 



To be more precise about these ideas, it is im- 
portant to set down their corresponding analytical 
expressions. The joint probability distribution for 
the Gaussian random processes x(t) and y(t) is, 

to/ n 7 7 dxdy Vx 2 +y 2 —2 P xy'] ( . 



where 
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When eq (2.1) is transformed to the polar coordi- 
nates (/?, <j>) of figure 1, and the phase averaged out, 
the amplitude distribution becomes: 
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This expression reduces to the Rayleigh distribution 
if the component correlation p is assumed to be zero. 
However, this assumption is probably not justified 
for electromagnetic scattering — as it is for electrical 
shot noise — and the more general expression (2.2) 
may explain interesting anomalies in the experimental 
envelope distribution data [10]. The data generally 
seem to emphasize larger values of R (relative to the 
mean) more than the Rayleigh distribution predicts, 
with corresponding reductions in the occurrence of 
small values. This behavior could also be explained 
by: (1) Admixture of a steady signal component — 
perhaps due to partial ionospheric reflections above 



the MUF, or (2) trends in the average signal, which 
would produce larger apparent deviations from the , 
mean level when viewed over a finite data sample. s 
On the other hand, such anomalies may represent 
the effect of meteor-reflected signals. 

One can also discuss the joint amplitude distribu- 
tion for two amplitude processes i?i and R 2 [5]. 
These could refer to signals received at two adjacent 
aerials, such as are often used in diversity receiving 
systems. Ri and R 2 could also refer to two time- 
displaced records of the same circuit signals, as 
shown in figure 1. Dual process distributions intro- 
duce additional statistical parameters /*, which 
represent the (time, space, etc.) correlations between 
the component processes Xi, yi, x 2 , y 2 . The purely 
probabilistic predictions of such distributions have 
been confirmed on tropospheric scatter circuits [11] 
but have yet to be studied carefully for ionospheric 
scatter. 

In general, one can say that the statistical predic- 
tions of scattering theory are consistent with the 
experimental data. This agreement is neither exact 
nor unique, but it is satisfactory. The real problem 
in scatter theory is to predict successfully the statisti- 
cal parameters, a, p, /x, etc., as functions of the 
variables which are experimentally controlled: radio 
frequency, path distance, time of day and season, etc. 
No theoretical estimates for the correlation coeffi- 
cients p and At have been published, with the excep- 
tion of certain calculations for tropospheric scatter 
[5], which may or may not be applicable to iono 
spheric circuits. Research to date has been focused 
exclusively on the average power, which according 
to eq (2.2) is 



<R 2 > = 2a 2 . 



(2.3) 



3. Electromagnetic Scattering by Turbulen 
Irregularities 

The basic wave equation which discribes radio 
wave propagation through an ionized plasma of elec 
tron density N with stochastic fluctuations 6N addec 
thereto is: 



[V 2 +k 2 -r (N +SN] E=0, 



(3.1] 



where k=2ir/\ is the free space wavenumber an(? 
r =2.8X10~ 13 cm the classical electron radius. Th; 
mean electron density N in (3.1) is usually omitted 
since refractive effects have been found to be negli 

gible [12]. In establishing expression (3.1), one as 

-> 
sumes that: (1) The stochastic variations 8N(r,t 
rearrange themselves in a time which is long com 
pared to the period of the radio oscillations, and (2i 
that gradients of 6N are negligible over a wavelength 
The latter assumption is clearly the more question 
able. 

It is not possible to solve eq (3.1) exactly, becaus 
6N(r,i) is an unknown stochastic function of position 
and time. Heretofore, the Born approximation ha 
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been used exclusively to describe single scattering by 
such fluctuations. This is essentially a three-dimen- 
sional iterative solution of the wave eq (3.1), which 
gives the scattered field strength at the receiver R 
in terms of the unperturbed wave E (r). 



E s (R)=r ( d¥ G (R,r) 6N (r,t) E (r), (3.2) 



where 



and 
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P T is the transmitted power and G T the transmitter 
gain function. The integration of (3.2) is extended 
over the common volume V illuminated by both the 
transmitter and receiver. 

Since both the transmitter and receiver are usually 
many wavelengths from the scattering blobs in the 
common volume, one can use the far-field approxi- 
mation [5] to express the scattered field strength as: 



E s {B)c 



°V4tT 47TJ/' 
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-6N(r,/) A / G'r (3.5) 



where the various quantities are defined in figure 2. 
The orthogonal random signal component x(t) and 
y(t) of figure 1 are identified with the real and imagi- 
nary parts of this expression respectively. The re- 
ceived power is computed from this expression by 
averaging its square and multiplying with the ef- 
fective receiver area. In terms of the scattering 
-^ > 

difference vector K=ki—k 2 '. 
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where A R is the actual receiver area. The space 
correlation probably vanishes unless r and r' are 
within several kilometers of one another, and the 
other factors are relatively weak functions of position, 
so that R 1 =R[, R 2 =R2, etc., to a good approxi- 
mation. Transforming the double integral to sum 
and difference coordinates (r,p), one finds: 



p «= p ^f r wiH K - f 



<*N(r f 0«N(r+p,*)> 

= PTrl jr^)/ r wi S{K) ' 



(3.7) 



where S(K) is the spectrum of turbulent irregular- 
ities to he discussed presently. This expression is 
probably a valid approximation, at least insofar as 
single scattering is the dominant scattering mecha- 
nism. 




Figure 2. Geometry for analysis of electromagnetic scattering 
from ionospheric irregularities. 

Transmission expressions which have been used 
thus far in interpreting experimental data acquired 
on forward scatter communication circuits, however, 
represent further assumptions, which arc 4 well worth 
identifying. 

a. The ionospheric irregularities are assumed to 
be isotropic, 3 so that the spectrum S(K) appearing 
in (3.7) depends only on the magnitude of the scat- 
tering difference vector. 



K=\K\- 



4tt . f0\ 

-Y 8in W 



(3.8) 



where 8 is the scattering angle between the upgoing 
and downgoing rays at the integration point r. The 
scattering process thus acts like a narrow band filter 
on the spectrum, emphasizing the uniquely important 
wavenumber given by eq (3.8). The isotropy as- 
sumption evidently precludes a discussion of mag- 
netically-alined irregularities with the resulting 
transmission expression. 

b. The volume integration (d 3 r) is collapsed by 
multiplying V by an average value of the integrand. 
The average scattering angle 6 is usually chosen to 
be that at the path midpoint, corresponding to the 
strongest scattering contribution Ji.e., smallest K). 
The effective volume is set by 6, the transmitter 
gain, and the width of an effective scattering la3 r er b 
(see fig. 2) : 



V-- 



---^- R\b esc-- 



(3.9) 



3 Anisotropic correlation functions were used by Staras in his extensive analysis 
of tropospheric scatter circuits and by Booker in discussing the backscattering 
from field-alined auroral irregularities. However, these techniques have not 
been applied to ionospheric forward scattering to date. 
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c. Polarization effects have been neglected, al- 
though a relatively simple theory predicts a factor 

-> 

sin 2 x in eq (3.6), where x is the angle between E 

-> 
and the scattering direction k 2 . More precise 
accounts [13, 14] of the vector field response support 
this approximation. 

The resulting transmission equation, 

P R =P T rlA R G R || esc (£) S(~ sin (3.10) 

has formed the only linkage between turbulence 
theory and experiments thus far. Expression (3.10) 
is evidently a substantial approximation to the actual 
electromagnetic response of the irregularities, and it 
is well to be aware of its limitations. On the other 
hand, a new technique has been developed [5] which 
allows one to remove all of the above approxima- 
tions — within the framework of single scattering — fit 
the expense of complicated integrations. 

The scattering cross section (per unit volume per 
unit solid angle) concept is often employed, viz, 



a(0,\)=T* o S 



/4tt . d\ 



(3.11) 



However, this approach does not emphasize the joint 
frequency-distance dependence of (3.10) which is 
characteristic of scatter propagation. In any case, 
it is quite clear that the spectrum of irregularities is 
the basic description of the turbulent electron density 
variations which is required by the electromagnetic 
theory. 

4. Turbulence Theories 

Results of the last section show that the spectrum 
of turbulent electron irregularities is the physical 
quantity of central interest for scatter propagation. 
More particularly, one is interested in the spectrum 
at wavelengths 



„ 4tt . 6 
K=— sin -; 

A. 2 



(4.1) 



which lie in the following range for typical vhf 
scatter paths. 



0.2m- 1 <K<1.3m- 1 , 



(4.1a) 



On the other hand, the wavenumber K s at which 
diffusion-viscosity effects ought to become impor- 
tant is within— or at least near — the range of 
interest. 

/ V 3 \ 1/4 
K.=(Pi) ^(Um- 1 (4.2) 






In terms of the spectral ranges identified in figure 3, 
this means that one is probably interested in the 
spectrum for the difficult transition region between 
the inertial subrange and the asymptotic dissipa- 
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Figure 3. Typical spectrum of dielectric irregularities plotted 
logarithmically versus wavenumber, identifying various 
turbulence ranges. 



tion range. This is unfortunate, since one probably 
cannot use the dimensional analysis methods which 
have proven so powerful in the inertial subrange. 
To avoid notational confusion, we use S(K) to de- 
note the three-dimensional spectrum, which is de- 
fined as the Fourier transform of the space correla- 
tion function. 



<«N(r,0«N(r+fi,0> 



"5? J 1 



d?ke ik ' R S(K). 



(4.3) 



If the turbulent irregularities are isotropic, the 

spectrum depends only on the magnitude of the 

-> 
vector wavenumber k } and the entire description is 
much simplified. 

The one-dimensional spectrum is denoted by T(Ic), 
and corresponds to wavenumber contributions to 
the mean square variation. 

<8N 2 >=^- 2 r dkk 2 S{K)= CdkT(k). 
^ Jo Jo 



This establishes the relation : 



(4.4) 
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Three distinct physical theories have been de- 
veloped in recent years to describe the irregularities 
of electron density established by the turbulent 
velocity field of neutral ionospheric gases: 

a. Theory of pressure fluctuations, 

b. Obukhov-Corrsin mixing theory, 

c. Villars-Weisskopf mixing theory. 

Each of these theories has now been extended into 
the dissipation regime, and the essential results are 
tabulated for convenience in figure 4, together with 
brief references to the many contributors. Without ^ 



Theorv 



Inertia! Range 



References 



Dissipation Range 



References 



1. Pressure *> -yl (15) Villars and / \4 (15) Villars and 

Fluctuations f o) k 4/3 fc -13/3 Weisskopf N o k 4/3 k 20/3 k -ll Weisskopf 

o V C o o \C J o s .... .,„ . 

v J K y (16) Wheelon 



dN 



2. Villars - 

Weisskopf 

Mixing V dh J 



(17) Villars and 

Weisskopf 

(16) Wheelon 
(19) Wheelon 

(18) Gallet 



'«oT 



k 8 k" 13 



Villars and 

Weisskopf 
(Unpublished) 

(16) Wheelon 

(19) Wheelon 



3. Obukhov- 
Corrsin 
Mixing 



dN 



k - 4/3 k- 11/3 



(20) Obukhov 

(21) Corrsin 

(22) Silverman 

(23) Bolgiano 

(24) Batchelor 



k" J D « v 

,-23/3 



D » v 



(24) Batchelor 

(25) Batchelor, 

Howells, 
and 
Towns end 



Figure -1. Summary of spectral models S(k) for electron density irregularities in the lower ionosphere. 



attempting to plead one case or another in the color- 
ful but unfortunate controversy which has been 
kindled in this new field, it is important to note 
several points which now seem rather clearcut: 

a. The theory of pressure fluctuations has essen- 
tially been abandoned as a major source of iono- 
spheric irregularities [17, 18]. 

b. The two mixing theories both agree that a 
(random) conveetive mixing-diffusion equation of 
the form: 



Q- r Dv^N(r,t) + 



v(r,tyvBN(r,t)=v-vN (4.5) 



characterizes the process, where v (r, t) is the turbulent 
solution of the Navier-Stokes equations. 

c. The Obukhov-Corrsin theory is patterned 
precisely on the Kologomoroff theory of the velocity 
spectrum. It essentially disregards the differential 
eq (4.5), assumes that the input and inertial range 
do not overlap, and arrives at the following result 
by dimensional analysis: 



T(k)^k-^^yh-^. 



(4.6) 



d. Villars and Weisskopf noted that the gradient 
of iV appears in the basic eq (4.5). Accepting 
dNc/dh as a legitimate parameter for dimensional 
analysis, one finds: 



r (&)-&- 



/</A r oY 
V dh ) 



(4.7) 



e. In deriving expressions for the transition and 
dissipation ranges, both mixing theories use approxi- 
mate analysis which are open to question. 

f. None of the theories discussed here have as yet 
been extended to include the influence of the earth's 
magnetic field. 



The reader is referred 
indicated in figure 4 for 
precise expressions. 



to the basic references 
detailed arguments and 



5. Signal Level and Scattering Heights 

One can combine the approximate transmission 
expression (3.10) with various turbulence models 
listed in figure 4 to predict theoretically the received 
power on a particular scatter link. This calculation 
and comparison has been performed for ionosphere 
forward scatter circuits so far [17, 26] only for the 
Villars-Weisskopf mixing model (4.7), which is 
subject to a rather stringent test since it contains no 
adjustable turbulence parameters. Using only the 
inertial subrange result for illustration, one has for 
the theoretical transmission equation: 



i\ 



bA B G B 



K£)'-(i)(f)" 



(5.1) 



All of the factors in this equation are known for a 
given circuit, except the appropriate value of 
electron density gradient. This evidently depends 
on the ambient electron density profile at the 
common volume altitude. 

In figure 5 we reproduce a suggested typical 
noontime profile of electron density for the lower 
ionosphere from reference [26]. The sharp lower 
boundary at 70 km is consistent with vlf reflection 
data, and the upper portions agree moderately well 
with rocket data. According to both mixing theories, 
large gradient portions of the profile ought to give 
rise to more intense turbulent irregularities and 
hence stronger scattering. A sharp gradient at 70 
km would reproduce the strong daytime signal 
contribution. The power levels measured on the 
Cedar Rapids-Sterling path [7] require that this 
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Figure 5. Possible noontime profile of electron density in the 
lower ionosphere which is consistent with various radio 
measurements from reference [26]. 



gradient be approximately 
-2,000 



dNo 
dh 



cm* 



per km, 



(5.2) 



which represents a very fast rise to the modest total 
electron density at this height. It would be interest- 
ing to see if the Obukhov-Corrsin theory gives a 
smaller number. 

A weaker gradient at approximately 85 km is 
suggested in figure 5, and should also contribute to 
the scatter signal. The gradient at this height is 
thought to be 



dNo 
dh 



-3C0 —z per km, 
cnr 



(5.3) 



corresponding to a signal some 10- to 15-db weaker 
than the 70-km contribution. It is interesting that 
both of these signal components are found in the 
short pulse measurements of Pineo [27]. In general, 
he finds unequal signals from two levels separated 
by about 15 km. Further evidence for weaker scat- 
tering from an upper layer com.es from the long range 
NBS experiments [28], which compared 36-Mc signals 
over 2270 km with 50-Mc signals received on a 
1240-km path. Earth screening on the longer path 
obscured all volume elements below 87 km in these 
experiments, and the (frequency-scaled) received 
signal was some 10-db weaker than the shorter 
circuit, which presumably could exploit both the 
70- and 85-km ranges. Further evidence for this 
dual stratification comes from the vertical incidence 
pulse backscatter measurements made at hf (see 
[26] for references). 

6. Diurnal and Seasonal Variations 

The signal level of a vhf scatter circuit exhibits 
large changes during each day, year, and sunspot 
cycle. Experimental data for the'NBS 50-Mc Cedar 
Rapids-Sterling circuit [7] is shown in figure 6, to- 



gether with the long range Newfoundland-Azores 
path. The first record is typical of most vhf paths 
and shows a 10-db rise from morning to noon, fol- 
lowed by a 15-db drop to postsunset values. The 
corresponding meteor count is strongly peaked at 
0600 and is essentially negligible between noon and 
midnight. This suggests that the afternoon and 
early evening is under strong solar control. 

A theoretical explanation [26] for the afternoon 
signal drop is offered by turbulent mixing theories 
through the mean profile gradient factors, since the 
mean profile is under very strong solar control all 
day. The noontime signal is contributed by the 
70-km level (sec. 5) and the afternoon decay shown 
in figure 6 ought to reflect the way in which this 
sharp gradient is erased by recombination at sunset. 
The rapid drop observed corresponds to the very fast 
recombination rates at lower levels. On the other 
hand, the decay of the upper gradient is found to be 
much more gradual, and in fact should produce 
measurable scatter signals all night. The noon-to- 
evening 1500-km signal drop in figure 6 presumably 
corresponds to the shift of scattering height from 
the 70- to 85-km levels, which is confirmed by pulse 
measurements [27]. The comparative absence of a 
steep signal drop on the long range path is again 
consistent with a slow recombination of the upper 
gradient. This line of reasoning forms quite a satis- 
factory explanation for the existing data on diurnal 
variations. The basic key to this agreement is the 
relationship between turbulent intensity and the 
mean profile provided by the mixing theories. 
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Figure G. Diurnal variation of power on two vhf scatter paths 
from reference [?]. 



Recent published data [29] on correlation of vhf 
scatter signal strength and sunspot number demon- 
strates a strong solar control over a 6-yr period. A 
good explanation is again provided by the relation- 
ship between turbulence intensity and the mean 
profile (gradient), since the latter is clearly under 
consistent solar control. 
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7. Frequency Dependence 

The basic transmission eq (3.10) depends upon 
the radiofrequency through: (1) The turbulence spec- 
trum S(K), and (2) the effective receiving area A R G R . 
To avoid confusion on how the second effect is ac- 
counted for, it is convenient to compare received 
power on two different frequencies for scaled anten- 
nas, viz, equal G R factors. The definitive experi- 
ments by the NBS group [10] on frequency depend- 
ence actually used carefully scaled rhombic antennas, 
so as to insure identical illumination of the scattering 
layers and similar antenna response. With this con- 
vention, the frequency variation is that due to the 
spectrum plus a X 2 factor which represents the physical 
scaling of the receiving aerial size. If the spectrum 
is expressible as a simple inverse power (n) of wave- 
number in the relevant turbulent range, one would 
predict a (scaled antenna) frequency scaling law [30], 



Pa-- 



const 



fn+2 



(7.1) 



The most recent experiment;! I data is that reported 
by Blair [10], using 30, 40, 50, 74, and 108 Mc. He 
finds that the exponent n in eq (7.1) is variable, 
ranging from 7 to 10. This variation is both diurnal 
and seasonal, so that it is dangerous to quote a 
mean value. During Late afternoon hours, however, 
when scattering from turbulence is presumably the 
dominant mode, the exponent n is definitely higher, 
and a value of 8 to 9 may be representative. Com- 
paring this with the turbulence results tabula led in 
figure 4, the Obukhov-Oorrsin inertial range result 
gives/ -56 , which is definitely too low. The Villars- 
Weisskopf inertial result f~ 7 is occasionally correct, 
but the corresponding dissipation behavior/" 15 is too 
high. The/ -96 dissipation result for d >> v given 
by Batchelor, Howells, and Townsend [25] describes 
some of the records. On the balance, however, one 
finds that none of the theories give a satisfactory 
explanation for the important frequency exponent 
variability. 

Earlier experiments of frequency dependence [7, 31] 
indicated a higher exponent scaling law between 50 
and lOOMc, than that measured between 30 and 
50Mc. Because the scattering process acts like a 
narrow band filter on the spectrum (sec. 3), it was 
suggested that the viscosity transition wavenumber 
of the ionosphere (4.2) did, in fact, divide the wave- 
number range of experimental interest (4.1). The 
steeper slope approaching the dissipation range and 
the smaller inertial slope seemed to explain the data 
rather nicely [16]. However, the newer data [10] 
indicates that such a dichotomy does not exist, and 
that the exponent is constant over the wavenumber 
range employed at any given hour — even though the 
slope itself shows diurnal and seasonal variations. 

8. Distance Dependence 

If the turbulent spectrum function S(K) is expres- 
sible as a simple inverse power (n) of wavenumber in 



the relevant range, two factors which depend on path 
length d occur in the basic transmission eq (3.10), 



Pr= 



constant 



/>l(sin|)" 



(8.1) 



If a denotes the mean radius of the earth, the distance 
R 2 from the average scattering point to the receiver 
is given by 

l?i=a 2 +(a+h) 2 -2a((i+h) cos CQ a/|Yfor a^>d. 

(8.2) 

The scattering angle is a more complicated function 
of distance d and scattering height h [16]. 



tan 



i—(iv 



""(£) 



2a) a d . 2k - ; ^ /0 , 
c^-r— —y for <7<Ca. (8.3) 
Aa d 



A graphical evaluation of B versus d for various h is 
given in reference [16]. By combining all of the 
above expressions, it is possible to establish the vari- 
ation of received power with distance predicted by 
any turbulence model (//). However, one should also 
note that earth screening horizon limitations have the 
effect of raising the lowest scattering midpoint as the 
distance is increased. Insofar as the strongest scat- 
tering centers are located at the lower height (sec. 6), 
this vertical inhomogeneity must be reflected in an ad- 
ditional distance dependence of the mean signal level. 
The experimental data now available is probably 
not accurate enough for a decisive test of the theories. 
The early data [7] used successive measurements of 
50-Mc signals at 491, 592, and 811 km, whereas one 
would have preferred simultaneous data over longer 
paths. The shipborne DSIR experiments [31] on 
41 Mc gave limited results in terms of defining the 
distance scaling law [30]. Recent (unpublished) air- 
borne receiver experiments made by the Lincoln 
Laboratory group out to 3000 km on 50 Ale should 
provide a valuable basis for testing the various theo- 
retical results. 

9. Sudden Ionospheric Disturbances 

A sudden increase in the flux of ionizing radiation 
from the sun increases the electron density rather 
quickly in the lower ionosphere, which in turn in- 
creases the nondeviative absorption of radio waves 
which pass through the D region [33]. This prompt 
absorption is called a sudden ionospheric disturbance, 
and affects both totally reflected short wave (hf) 
signals and the partially reflected waves responsible 
for (vhf) ionospheric scatter propagation. A typical 
set of simultaneous signal amplitude records made on 
a shortwave circuit at 6 Mc and two adjacent scatter 
links at 28 and 50 Mc is reproduced in figure 7, from 
reference [7]. The shortwave link drops out im- 
mediately, while the 50-Mc scatter link increases 
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Figure 7. Transient behavior of power on three frequencies 
taken from reference [7], 



gradually until the recovery cycle of the transient D 
region begins. The 28-Mc link shows a mixture of 
the two effects. Viewed in a qualitative way, these 
data lend strong support to the notion that iono- 
sphereic scatter links are under direct solar control 
much of the time [29]. 

One can construct 4 an analytical explanation for 
the transient amplitude behavior shown in figure 7. 
Nondeviative absorption of a radio wave of frequency 
f can be represented approximately bv a multiplica- 
tive factor [33] ; 

exp-^j^dyN(y,t)v(y) } (9.1) 

where h is the reflection height. N (y,t) is the in- 
stantaneous value of electron density at the ray path 
integration point y, and v is the collision frequency 
at the same point. The totally reflected shortwave 
field strength is modified by the factor (9.1) alone, 
which is very large because the carrier frequency / 
is relatively small. The rapid realization of this 
absorption loss indicated by figure 7 indicates how 
rapidly the electron density profile is modified by the 
burst of solar ionization. 

If one is to believe the turbulent mixing theories, 
the field strength of an ionospheric scatter circuit 
depends both on the transient electron density profile 
and its gradient, viz, 



P~ 



r dN(h, t) 
L dh 



Jexp- l£dyN(y, t)v{y). (9.2) 



* A. D. Wheelon, unpublished calculations. 
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Because of height-dependent recombination rates, 
etc., rapid changes in the flux of ionizing radiation 
produce substantial modification of the basic elec- \ 
tron density height profile. Careful analysis shows 
that transient modification of ionization gradients 
in the D region can explain the signal increase of 50 
Mc shown in figure 7. (The exponential factor is 
negligible at such high frequencies.) On the other 
hand, the 28-Mc signal ought to suffer significant 
absorption, and the competition between enhanced 
gradient scattering and this absorption can be ex- 
plained with the theory. A possible explanation for 
this type of behavior is thus provided by turbulent 
mixing theories and standard absorption expressions. 
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